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Different mathematical methods have been applied to obtain the analytic result for the massless 
triangle Feynman diagram yielding a sum of four linearly independent hypergeometric functions of 
two variables F±. These are defined for especific regions of convergence for the ratios of the squares 
of momentum variables. In this paper I work out the diagram and show that that result, though 
mathematically sound, is not physically acceptable when it is embedded in higher loops - meaning 
further momentum integrations - because it misses a fundamental physical constraint imposed by 
the conservation of momentum, which should reduce by one the total number of linearly independent 
(Li.) functions F4 in the overall solution. Taking into account that the momenta flowing along the 
three legs of the diagram are constrained by momentum conservation, the number of overall Li. 
functions that enter the most general solution must reduce accordingly. 

To determine the exact structure and content of the analytic solution for the three-point function 
that can be embedded in higher loops, I use the analogy that exists between Feynman diagrams and 
electric circuit networks, in which the electric current flowing in the network plays the role of the 
momentum flowing in the lines of a Feynman diagram. This analogy is employed to define exactly 
which three out of the four hypergeometric functions are relevant to the analytic solution for the 
Feynman diagram. The analogy is built based on the equivalence between electric resistance circuit 
networks of type "Y" and "Delta" in which flows a conserved current. The equivalence is established 
via the theorem of minimum energy dissipation within circuits having these structures. 

PACS numbers: 11.10.Gh,11.15.Bt,11.55.Bq, 12.38.Bx, 14.70. Dj 

I. INTRODUCTION 

In working out the higher order loops of Fcyman diagrams, one has to perform a growing number of momentum 
integrals. It may happen that the total integration may be performed via sectioning it so that first "simpler" lower- 
loop integrals are performed and then the result embedded in a higher order diagram integral. This method of doing 
things is not new and have been applied ever since higher loop calculations have been performed in quantum field 
theory calculations. There are several types of one-loop integrals that occur in these more complex integrands of 
higher loop diagrams; among which we cite the one-loop self-energy type integrals and one-loop triangle diagram type 
integrals. In this work I focus my attention in the latter type of integrals. 




FIG. 1: One-loop triangle Feynman diagram. 



Massless triangle Feynman diagrams (see FigJT]) occur in field theoretical perturbative calculations, being one of 
the primary divergent (one-particle irreducible) graphs for the non-Abelian gauge theories such as Yang-Mills fields. 
This makes this kind of perturbative calculation the more interesting and necessary once we want to check the higher 
order quantum corrections for many physical processes of interest. 

Two decades ago, Boos and Davydychev [l[ were the first to obtain the analytic formulae for the one-loop massless 
vertex diagram using the Mellin-Barnes complex contour integral representation for the propagators. Later on, 
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Suzuki, Santos and Schmidt reproduced the same result making use of the negative dimensional integration 
method (NDIM) technique of Halliday and Ricotta Q for Feynman integrals, where propagators are initially taken to 
be finite polynomials in the integrand, and once the result is obtained, analytically continued to the realm of positive 
dimensions. 

These two widely different methods of calculation yielding the same answer lends to the analytic result obtained 
a degree of certainty as far as the mathematical soundness is concerned. The dimensionally regularized {D = 2u>) 
massless triangle integral reads 

I A (p 2 , q 2 , r 2 ) = J (fc2+ze )p_ p) 2 + ^p._ 9) 2 +le ] (!) 

in which p and oare two of the independent external momenta, and r = q — p the momentum for the diagram's third 
leg. It yields 0,12] 

Ia = tt" ( P 2 f ~ 3 {AF 4 (a, p, 7 , 7'; x, y) (2) 

+B X 1 -* F 4 (a + 1 - 7 , p + 1 - 7) 2 - 7 , 7 '; x, y) (3) 

+Cy 1 -'*' F±(a + 1 - 7 ', p + 1 - 7 ', 7l 2 - 7'; x, y) (4) 

+Dx 1 -V" 7 ' F 4 (a + 2 -7-7', £ + 2-7-7', 2 - 7, 2 -7'; x, y)} (5) 
where we have introduced the following definitions for the coefficients: 

r 2 (^-2)r(3-^) 
r(2w-3) 
B = c r(2-^)r(c-2)r(c-i) 
r(2w - 3) 

d = r 2 (2-w)r(cj- 1), (6) 

and the various hypcrgcomctric function parameters and variables are given by: 

a = 1, 

/3 = 7 = 7' = 3- w, 
^2 

2/ = %■ (7) 

As it stands by itself, the analytic result ©-(O is perfectly sound and good. However, when we consider it as 
a part of a higher order diagram calculation, it is no longer valid since the variables in the F^s are constrained by 
convergence rules which further momentum integration violates. So we need to look for the solution that will allow 
such further integration without violating convergence constraints. 

To begin our search let us consider the structure of the solution as it stands with four F^s. From the study of 
higher transcendental functions, consider the following system of coupled, simultaneous partial differential equations 

i 

,d 2 Z 2 d 2 Z n d 2 Z , „ „, .dZ , n dZ nr7 

x{\ - x) — - y 2 — - 2xy—— + [7 - (a + (3 + l)x} — -(a + p + l)y— - a(3Z = 
ox z oy z oxoy ox ay 

pp 7 pp y ppy p\y pp7 

- - x ~q^ - 2x y^o- y + ft' -(«+/* + -(«+£+ i)^ - »pz = o, (8) 

with variables x and y independent. The most general solution for such a system of differential equations consists of 
a linear combination of four hypergeometric functions of type -F4, as follows Q: 

Z = AF 4 (a, 0, 7, 7'; x, y) (9) 

+ Bx 1 -T^ 4 (a + 1 -7, £ + 1-7, 2-7, l'\x, y) (10) 

+ Cy 1 -T> 4 (a + l-7', p + l-Y, 7, 2-7'; x, y) (11) 

+ Dx 1 -^y 1 -T'F 4 (a + 2 -7-7', £ + 2-7-7', 2-7, 2-7'; x, y). (12) 
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Naive comparison between (J2j — (JSJ an d (j9|)- ([T2l) would tempt us to write Z = without much a do. However, 
we have to remember that for the case /a, variables x and y are not independent ones, since they are connected 
by momentum conservation: r = q — p. Therefore, the most general solution for the underlying system of coupled 
equations that are associated to the triangle diagram should not have four independent functions F$ but only three, 
on account of the one constraint due to momentum conservation connecting variables x and y. 

From a pure mathematical point of view then, the solution for the triangle Feynman diagram must obey this 
requisite, since not all the four hypcrgcomctric functions in it are linearly independent to each other. 

Therefore, on a sound mathematical basis, when we are interested in processes requiring momentum integration, this 
constraint must be considered and the result for the triangle diagram integral should be expressed as a combination 
of only three linearly independent hypergeometric functions F4. In other words, the result for triangle diagrams that 
is to be embedded in higher order diagrams must allow for all the momentum integration range, which is forbidden 
in the four term F4S format, due to the variables convergence constraint. 

Neither the Mellin-Barnes formulation of Boos and Davydychev nor the NDIM technique employed by Suzuki et 
al to evaluate the massless triangle graph take into account the basic physical constraint imposed on such a diagram 
due to the momentum conservation flowing in the three legs. Momentum (that is, energy and three-momentum) 
conservation is one of the basic, fundamental tenets of modern natural sciences, and plays a key role in determining 
the final form for the analytic result for the massles one-loop vertex correction without convergence constraint forbiding 
further momentum integration. 

One of the simplest examples where the four-term hypergeometric combination @ ~© above mentioned does not 
reproduce the correct result is when it is embedded in a higher two- loop order calculation Q. Consider, for example, 
the diagrams of Fig [2] 




FIG. 2: Two-point two-loop Feynman diagrams: Fying-saucer side-view and front-view. 



These two two-loop Feynman diagrams have similar loop momentum integrand structures, namely, 

d D kd D q 



r(2) _ 

L FS ~ 



(k 2 ) a q 2 [(k-p) 2 ][(q-ky 



(13) 



in which we omitted the infinitesimals (ie) associated to each propagator in the denominator and a can be either 1 or 
2 depending on which diagram one considers. 

We can perform the double integration in two ways. The easiest one is doing the q integration first and then the k 
integration, 



r(2) _ 

1 FS ~ 



d D k 



d D q 



{k 2 ) a [(k - p) 2 } J q 2 [(q-k) 2 } 
so that both integrands are of the one-loop self-energy type, and the results are: 

^2c-4 r3 (^ - !) r ( 2 - - 3)r(4 - 2u>) 



j(2) _ 2u 
1 FS — n 



[P 



(P 



r(2w - 2)r(3 - w)r(3w - 4) 

2u _ 5 r 3 (cj - l)r(2 - ^)r(2u; - 4)r(5 - 2lj) 



for a — 1, 



for a 



r(2w - 2)r(4 - w)r(3w - 5) 

The result (fT"5)) should also be obtainable from doing k integration first and then q integration at last: 



r(2)_ fdPq f 

FS ~Jq 2 J {k 2 Y 



d D k 



[{k-p) 2 ][{q-kf 



(14) 



(15) 



(16) 
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It is clear that this path is harder to travel, since the first integrand is of the one-loop triangle diagram type, (QJ. 

It is here that the plot thickens, because those results in (fT5"]) can never be achieved by using the triangle diagram 
as expressed in terms of (J2]) - ([5]). For neither of the diagrams the four-term solution for the triangle when embedded 
into the two-loop graph yields the correct answer. The right result can only be achieved when we use a three-term 
triangle solution embedded into the two-loop structure [|[ . What is the problem? I argue that the problem stems from 
the fact that these four terms ([2]) - (JSJ are not correct because they are not linearly independent solutions; physical 
momentum-energy conservation law must be considered, and this must reduce by one the overall number of F 4 's that 
should enter into the most general solution for Ia in equation (JT|). 

How do we go about in reducing those four terms into three? We could use the analytic continuation formula for 
the F 4 function 

/ , s T(d)T(b-a), , , z 1, 

F 4 (a,b,c,c';z,w) = ) ' ) '-{-w) a F 4 (a, a + 1 - d, c, a + 1 - b; -, -) 

I (o)l (a — a) w w 

+ S " M (^)'^(fe + 1 ~ d, b, c, b + 1 - a; £ , 1) (17) 
I (ajl (a — b) ww 

in order to reduce the number of independent functions from four to three. 

However, for the overall four hypergeometric functions F4, that appear in I a, one has at least two possibilities to 
pair them two by two according to (|17l) : Either combining and (@| or combining and ([5]). 

Making a suitable change of variables and redefining conveniently the various parameters of the hypergeometric 
function, we have, for example, pairing (j3|) and ((5|) according to (|17|) : 

x 1 -T{BF4(a + l-7,/3 + l-7,2-7, 7 '; a; ,y) + Dy 1 ^'^ + 2 - 7 - 7 ', /3 + 2 - 7 - 7 ', 2 - 7, 2 - 7 '; y)} 



fi(B, D)^ 4 ( P + 2 - 7 - 7', 1 + P ~ 7, 2 - 7, 1 - a + /3; -, - J 

V y 2// 



(18) 



where f2(B, D) is a new coefficient which depends on the previous ones B and D. 

Another possibility would be to make a similar suitable change of variables and redefinition of parameters in a 
convenient way so as to pair ([2]) and (|4]) according to ([T7| . Either way we accomplish the reduction to three linearly 
independent F 4 s without convergence constraint that forbids further momentum integration, but which one must we 
choose? 

In other words, which three of those four F±s contained in the result 2a should be considered cannot be determined 
by momentum conservation only; it needs another physical input to this end. 

In order to determine the exact structure and content of the analytic result for the triangle diagram I employ the 
analogy that exists between Feynman diagrams and electric circuit networks. The conclusion of the matter can be 
summarized as follows: since there is a constraint (call it either a initial condition or a boundary condition) to the 
momenta flowing in the legs of a triangle diagram, it means that the overall analytic answer obtained via Mellin-Barnes 
or NDIM technique, which comes as a sum of four linearly independent hypergeometric functions of type -F4, in fact 
must contain only three linearly independent ones. Which three of these should be is determined by the momentum 
conservation flowing through the legs of the diagram, and, by the circuit analogy, by the equivalence between the 
"Y-type" and "A- type" resistance networks through which conserved electric current flows. 



II. ELECTRIC CIRCUITRY AND FEYNMAN DIAGRAMS 



Although known for a long time and mentioned sometimes in the literature of field theory, the analogy between 
electric circuit networks and Feynman diagrams, more often than not stays as a mere curiosity or at a diagrammatic 
level in which the drawing is only a representation that helps us out in the "visualization" of a given physical process. 

There is, however, a few exceptions to this. The earliest one that I know of is by Mathews Q as early as 1958 
where he deals with singularities of Green's functions and another one by Wu in 1961, where he gave an algebraic 
proof for several properties of normal thresholds (singularities of the scattering amplitudes) in perturbation theory. 

What I am going to do here is just to state the well-known result for the equivalence between "Y"- and "A-circuits" 
which is used, for example, in analysing the "Kelvin bridge" in relation to the more commom "Wheatstone brigde" . 
The equivalence I consider is for the network of resistors, where energy is dissipated via Joule effect. Then with this in 
hands I argue which three out of the overall four hypergeometric functions F4 present in the one-loop vertex obtained 
via Mellin-Barnes or via NDIM should be considered as physically meaningful. 

Note that the equivalence between those circuits of resistors can be established based on the conservation of energy 
flowing in their internal lines. 
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FIG. 3: "Y-type" and "Delta-type" resistor networks. 



Consider then the "Y"- and "A-circuits" of Fig H 

These two network of resistors are related to each other; for example, the "Delta-circuit" can be made equivalent 
to the "Y-circuit" when resistors obey the following equations [|| : 



Ra 
Rb 

R c 



R2 

Ri 
Ri 



Rs 
R3 
R 2 



R2R3 

Ri 
R1R3 

R~2 
R1R2 
R~ 3 



(19) 



In electric circuits a general theorem guarantees that the current density in a conductor distributes itself in such a 
way that the generation of heat is a minimum Q. Strictly speaking, since I consider the "Y"- and "Delta- networks" 
by themselves and not embedded within an electrical circuit, I cannot consider the minimum of power generation, but 
still I can consider the equivalence between the power generated in one of them — say "Y-circuit" — as compared to 
the power generated in the other one, the "Delta-circuit" . 

Labelling the currents that flow through resistors R2 and R 3 of the "Y-circuit" by p, q and r, respectively, by 
reason of current conservation we have that, for example fioj 



P 



(20) 



Each of these currents will generate heat according to Joule's law and the corresponding power will be given by 
Ohm's law. So in each leg, the product of the resistance times the square of the current will give the power generated 
by the current flowing through it. The overall power generated in the "Y-network" is therefore: 



P 2 Ri 



q 2 R2 



r 2 R 3 



Ptotal 



(21) 



Dividing (|2"Tj) respectively by p 2 , q 2 and r 2 we have the following system of equations: 



Ri 4 

i*. 

V 
X 



yRil 

I- -i?2 

Z 



xR 3 
zR 3 
+ R3 



n -^total 
= -Ptotal 

q* 



— ^--Ptotal 



— R r 



— R n 



(22) 
(23) 
(24) 



where for convenience I have defined z = — = — and R p = 



Pi 



total ■ 



Multiplying the last equation ([24]) by x and comparing with (|22|). 



Ri 

Ri- 



-yR2 
.r 



XR3 — —5 ^total 



-R 2 + xR 3 



total 



Rn 



R, 



(25) 
(26) 
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x 

it follows that the coefficient of R 2 is y = — . 

z ^ 

In a similar manner, multiplying (|22[) by — and comparing with (|23[) 



I x 1.1 

—Ri + R 2 H — i?3 = — ^-Ptotai = —R-p (27) 

2/ y yp y 

-Ri + R 2 + zR 3 = -iPtotai = -R P (28) 

y t y 

X 

it follows that the coefficient of R 3 is z = — . 

V 

Finally, multiplying (|24[) by yz and comparing with (|22[) , 

R 1 +yR 2 + xR 3 = ^P tota i=R p (29) 

— i2i + y-R 2 + = ^Ptotai = — R P (30) 

it follows that the coefficient of R\ is 1 = — , showing the thorough consistency among these results. 

x 

On the other hand, in the "Delta- network" , let I a label the current flowing through resistance R a , so that the power 
generated in this resistor is, using (fT9|) 

P a = I 2 R a = I 2 + ^2 + ^3} (31) 

Now, as the resistances follow the equivalence defined by ((TO)) , the current I a that flows through R a must be such 
that: 

I 2 a R a = p 2 a R ha + ql R 2 , a + r 2 a i? 3 ,a (32) 

where 

R 2 R 3 

Rl - a ~ ~rT 

R2,a = R2 (33) 
R3,a = R3 

so that the corresponding currents flowing through them must be "scaled" accordingly as 

pI = v -p 2 =p 2 

x 

<& = yi 2 (34) 

r 2 a = zr 2 (35) 

Therefore, the solution for the one-loop triangle Feynman diagram must be such that it contains three terms which 
are proportional to "currents" p 2 , yq 2 and zr 2 . This solution is exactly achieved when I combine © and (O in a 
similar way as was done in (|T8)) . Combination of ([2]) and (J4j) as in (|T8|) leaves the term ([5]) in the solution, which is 
proportional to xyp 2 , and therefore not suitable since it violates the equivalence above demonstrated. 

Explicit analytic solution for the one-loop massless triangle Feynman diagram reads therefore 

£ nc-ioo P (p, <?' r ) = 7r"(p 2 r" 3 r p .F4(l,3- - oj,3 - oj;x,y) 

+ 7T"y(q 2 r- 3 T q F 4 (l, w - 1, 3 - u, oj - 1; x, y) (36) 

/ x V 

+ 7r"z(r 2 )"- 3 r r P 4 1, w- I.lj- 1,3 -lj;z= -, - 

V y y 

where 

r __ r __ r _ r 2 (^-2)r(3-^) 

p_ q_ r r(2w-3) 

with a; being the dimensional regularization parameter. One notes that the variables for the last -F4, though related 
to those that appear in the first and second F48, differ from them. 

This construction is particularly important for embedded triangles of higher order Feynman diagrams, since the 
expression (|36p bears in it two- variable hypergcomctric functions of different types, namely, x with y and z with y^ 1 
which covers different regions of convergence. 
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III. CONCLUSION 



We have seen that the previous result given in the literature for the massless one-loop vertex correction of Feynman 
diagrams in a general field theory /a with four linearly independent hypcrgcomctric functions of type F± can be 
compared to the most general solution of a system of coupled differential equations for higher transcendental functions. 
These pressuposes that the variables are independent ones. Since the result Ia has exactly the same structure as the 
most general solution for the system of differential equations, it follows that the result for the vertex correction can 
be thought of as having a underlying system of differential equations that it must obey. However, contrary to purely 
mathematical system of differential equations of two variables, the momentum flowing through the Feynman diagrams 
obeys a physically important conservation law, which connects the variables x, y and z among themselves, via, e.g., 
r = q — p relation. So, the most general result for it should not have four independent two variable hypergeometric 
functions F4, but only three. 

Using the basic physical principle of momentum conservation I deduced the correct analytic result for the one-loop 
massless triangle Feynman diagram in terms of three linearly independent hypergeometric functions F4. The reason 
why there must be only three linearly independent functions in the solution is due to the fact that the variables x and 
y are such that implicit in them is a momentum conservation constraint that must be taken into account properly. 
The set of which three linearly independent functions F4 should be is given by another physical input, deduced in 
analogy to the equivalence between "Y"- and "Delta- network" of electrical resistance circuits. With the conservation 
of momentum correctly taken into account, the final result for the mentioned Feynman diagram is the physically 
correct and relevant analytic solution to the problem, having only a combination of three linearly independent F4 
functions in it. 

Moreover, of particular importance is this three-term result when vertex corrections are inserted in more complex 
Feynman diagrams. If one takes the four-term result I a with all the hypergeometric functions with same variables 
x and y, when this is further integrated in loop momentum k of higher order diagrams, it cannot be integrated 
on the whole domain of integration — 00 < k < +00 since all the F^s entering the solution /a has the same and 
definite domain of convergence. Not so with the (|36[) where the third F4 allows for the complementary domain region 
of momentum integration. To understand spcciffically which regions of convergence are for each of the two-variable 
hypcrgcomctric functions F4 see for example [J]. Having said that, whenever one has to deal with calculating multiloop 
Feynman diagrams in which the one-loop vertex correction is a subdiagram thereof, the correct result to embed in 
the integrand of the multiloop diagram integral must be the one given in (|36j) . 

Finally, for the sake of completeness, I would like to emphasize that two other symmetric results can be obtained just 
by cyclic permutations among (p, q, r) momentum variables, following the original symmetry of the triangle diagram. 
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